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ABSTRACT 


An experimental procedure is presented for lifetime testing of 
graphite bundles under constant load. The attributes of the experiment are 
expedience in implementation and a_ substantial accumulation § of 
information equivalent to a large number of single filament tests. To © 
achieve the objectives of the experiment a specially adapted Instron 
machine was used with a digital process/controller. Two trial tests were 
conducted using Hercules high strength graphite. The preliminary results 


are presented and the effects of inter-fiber friction evaluated. 
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I. INTRODUCTION 


Composite materials have been used by man for a very long.time. The 
first to be used were naturally occurring composites, such as wood. Then 
man discovered that there were advantages to be gained by combining 
materials with one component being fibrous, such as the use of straw to 
strengthen mud bricks. More recently, fiber reinforced composites with 
polymer and metal binder that have high strength to weight and high 
stiffness to weight ratios have become important in applications such as 
aircraft and space vehicles. 

In the 1960’s the U.S. Air Force began programs to explore aircraft 
structures made of composites. The first flight-worthy component 
produced was the horizontal stabilizer for the F-Il. Another major 
milestone was the production of a composite stabilizer for the F-14. That 
was followed by the composite stabilator for the F-15, and a composite 
rudder and stabilizer for the F-16. [Ref. 1] 

The advantages of composites when compared to conventional 
materials is that they exhibit the best qualities of their constituents and 
sometimes qualities that neither constituent possesses. Therefore, an 
understanding of reinforcement processes together with failure 
mechanisms is important for the development and production of high 
quality composites and for forecasting the long term stability of such 


structures, particularly in high performance applications. 


The purpose of this thesis is to present an experimental procedure 
for lifetime testing of graphite bundles under constant load for use in 
forecasting the life of graphite composite structures. And to demonstrate 
by comparison to single filament data that composite life is bracketed by 
composite load-sharing life as the upper bound, and fiber life as the lower 


bound. 


ll. BACKGROUND 


A. FIBER STRENGTH 

Brittle fibers used as reinforcement materials for advanced 
composites, such as carbon and qlass, are characterized by high strength and 
high modulus, and have tensile strengths that are statistical rather than | 
deterministic in nature. This statistical dispersion in strength is attributed 
to the existence of fiber surface defects. 

Rosen [Ref. 2 and 3] presented the results of an analytical and . 
experimental composite failure study. His theory considers fibers having a 
statistical distribution of flaws that result in individual fiber breaks at 
various stress levels. The load ina broken fiber is assumed to be distributed 
equally among the remaining unbroken fibers in a cross section. Composite 
failure is hypothesized to occur when the weakest cross section is unable to 
sustain the applied joad. 

In Rosen’s model, the applied load is considered to be supported 
entirely by the fibers since their extensional modulus is much greater than 
that of the matrix. Because the fibers have randomly distributed flaws they 
break randomly throughout the length of the fiber as the load is increased. In 
the vicinity of a break, the fiber stress is redistributed by the matrix binder, 
limiting only a portion of the fiber not fully effective in resisting the 
applied load. By a logical extension of this reasoning, the composite can be 
thought of as consisting of a series of layers of elements, analogous to links 


In a chain, whose axial length is equal to the ineffective length 4. 


2 


The ineffective length which is a function of the fiber volume 


fraction can be expressed, based upon an elastic analysis, as 


pt a 


8 =) ove? - NEt/G ql “cost {ll + (1+ 9/201 - Pde (1) 


where Vr is the fiber volume fraction, E, is the fiber Young’s mudulus, G,, is 


the matrix shear modulus, 9 is the fraction of the undisturbed stress value 


below which the fiber is considered to be ineffective, and dr is the fiber 


diameter [Ref. 4]. 

Although Rosen’s model agrees qualitatively with experimental 
data, in that failure is associated with the accumulation of fiber breaks, 
there is a disparity between predicted and observed failure loads, due to the 
simplifying assumption that the load in a broken fiber is uniformly 
distributed among the other fibers in the layer. 

Since it was logical to study the correlation between the theoretical 
strength of the weakest fiber and the observed failure loads, Rosen 
considered a population of fibers of length L whose strength is characterized 
by the probability density function g(o) and cumulative distribution function 
G(o). (Ref. 2 pp 68-71] 

For a sample of n links, where n = L/8, the probability density 


function for the strength of the weakest fiber is given by 


(0) = g(o)/n Itt - G(a)6/)-1 (2) 


It is assumed that the fiber strength can be characterized by a two 
parameter Weibull distribution of the form 


G(o) = I- exp(-«Lo8) (3) 


which has a corresponding density 


g(o) = xl BoP"! oxp(-xLo¥) (4) 


With L = §, and substituting equation (3) and (4) into equation (2) and then 
differentiating, the desired expression is obtained for the statistical mode 


of the weakest fiber strength distribution which is 


0” = Vi(xcb Bey 1/8 (5) 


The constants « and § denote the scale parameter and shape 
parameter of the Weibull distribution, respectively, and can be evaluated by 


using experimental strength-length data. 


B. FBeRe aire tine 
Consider a bundle of stiff, brittle fibers impregnated with a flexible 
matrix. Suppose a constant, tensile load is applied such that the impregnated 


bundle, though surviving at first, fails after many hours. This type of- 


fracture is a thermally activated process, and is referred to as 
stress-rupture or creeprupture. As athermally activated process, failure 
originates in the fibers where the molecules undergo random, thermal 
vibrations in time. As molecules slip or rupture, neighboring molecules 
become overloaded, thus increasing their failure rates. Such molecular 
failure accumulate locally and give rise to growing microcracks. These 
cracks eventually lead to broken fibers scattered throughout the composite. 
The randomness of these fiber breaks, both in position and time, is magnified 
by randomly distributed structural imperfections. [Ref. 5 pp 135-136] 

Historically, various kinetic modeis have been proposed to explain 
the phenomena of creep-rupture. But none of those proposals dealt with the 
variability of the lifetime data, which is an important aspect of 
creep-rupture in composites. 

Wagner (Ref. 6] presents the results of an analytical and 
experimental study for creep-rupture lifetime of Keviar 49, fibers. The 
model, reflects fiber variability and sensitivity of lifetime to load level, 


based upon the logarithmic approximation 


U(ad) = -Ug In(o/do) (6) 


where U(co) the thermal activation energy to rupture an atomic bond is a 


function of molecular stress o, the maximum bond force oo, and the 


activation energy Uo in the absence of stress. This approximating function 


|? 


i 


accounts for the power law relationship for the dependence of lifetime on 
stress level [Ref. 5). 

The power law framework which is based on the thermal activation 
process of molecular failure generates the result that the fiber strength and 
fiber lifetime follow the Weibull distribution. 


The Weibull distribution for fiber lifetime under constant stress is 


F(t) = 1 - exp{-[t/t)(o)]9} (7) 


with shape parameter 8 and scale parameter 
tho) =m YF te(0) (8) 
where m= 1/8 ane 
tsa) = o WUo/(KTIx,.1/8 (9) 
% iS apositive constant from the power law relationship. 


C. COMBOsiie LIFETIME 
As presented earlier, Rosen assumed the load in a broken fiber to be 
uniformly distributed among the other fibers in the layer which are intact. 


However, a more complex structure exists in the relationship between the 


fiber and matrix. Which is, that the fibers adjacent to a broken one are 
subjected to a load intensity greater than that which is sustained by 
fibers distant from the fracture point. Phoenix and Wu [Ref.5] discribe this 
phenomena as local load sharing. 

In local load sharing, the matrix serves two important functions: 
First, the effect of the fiber break is isolated longitudinally along the fiber 
as the shear stress in the matrix allows the fiber stress to return to normal 
a short axial distance away from the break. At the same time it permits the 
lateral transfer of the failed fiber load to its nearest neighbors. The 
increased loads on these neighbors greatly enhance their rate of failure. As 
these neighboring fibers break a transfer of the load continues until some 
unstable group of adjacent breaks, called a k’-crack emerges, and suddenly 
propagates across the composite. [Ref. 5] 

At high load levels, the shear in the matrix at fiber fractures may 
exceed the matrix mechanical properties causing cracks which propagate 
longitudinally as well as transversely. The ability of the composite to 
support high tensile loads is strongly influenced by the shear carrying ability 
of the matrix and by the ability of the fiber to sustain high tensile loads. 

The mathematical model of the failure process presented by Phoenix 
and Wu Is a weakest-link arrangement of independent bundles of length equal 
to the effective load transfer length &§ which is given by equation (1). 
In any given bundle the fiber elements share the applied load, which yeilds a 
nominal fiber stress. Since a failed element supports no stress, the stress 
of the failed element is redistributed onto its nearest neighbors, one on each 


side. An intact element next to one or more consecutive broken elements 


Carries a stress KX which is expressed as 


Ke=MM(2) +22) +, jz. (10) 


where K, is called the load concentration factor, X is the nominal fiber 


stress, andr is the number of consecutive broken fiber elements. 
The distribution for composite lifetime has the Weibull 


approximation 


Hm nt) * 1 - expl-[t/t (Lyk 9}, t > 0 (11) 


with shape parameter k*8 and scale parameter 


tel) = (mn) /K" 8) ayy te(L) (12) . 


where dyx is given by 


dx = T(B + 1-178 (kB + 1) (KB) Al - k*)/(k*B) MK; ) p/k* (13) 


The Weibull approximation given by equation (11), is valid for all values of 


Bp greater than 6, [Ref. 5 pp 141 - 142]. 


In these expressions the Weibull shape parameter and _ critical 
crack size varies with stress level, and account for the variability and 


size effect in strength and lifetime. 


Il. EAPERIMENTAL METHOD 


A. MATERIALS AND PREPARATION 
The fibrous material used to validate the experimental procedure to 
be presented was Hercules Magnamite high strength graphite with the 


following specifications: 


Type AS4 with W sizing 
Tow 3000 filaments 
Denier 0.005746 grams/inch 


Diameter 7.0 microns 


TO prepare the samples for the experiment, 100 inches of fiber was 
layed out on a table. One end of the fiber was secured to the table while 
approximately 5 inches of the other end hung over the table edge with a 
2 kilogram weight attached. A pulley at the edge, that the fiber layed over 
elimenated friction between the fiber and the table edqe. 

The 2 kilogram weight (approximately equal to 10% of the tow 
breaking load) was used to provide a small amount of tension and to ensure 
that all the filaments in each sample would be the same length when cut. 

A smear, about an inch in length of Tru Bond 5 Minute Epoxy Gel was 
applied to the fiber every IS inches. When the epoxy had cured the fiber 
was cut to produce six samples, each 15 inches in length, with one end of 
the sample coated with epoxy gel bonding its filaments together. Each 
sample was placed into a glass funnel that had a 10 inch exit tube, 


that protected the sample from damage during handling. Each funnel 
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with its sample was placed into a bath of methyl ethyl ketone for a period 
of twenty-four hours to remove the sizing, and then placed into an ethyl 
alcohol rinse (Figure 3.1) for one hour to flush away the residual sizing, 
and reduce random slack and disentangle the fibers. To prevent swirlling 
of the alcohol through the funnel during the rinse, the copper tube (Figure 
3.2) which delivered the alcohol into the funnel’s mouth, was designed with 
perferations about its perifery. Immediately after the rinse while the 
sample was still wet, to take advantage of the surface tension of the fluid, 
holding the fibers together and straight, the free end of the sample was 
smeared with the epoxy gel. The sample remained in the funnel until the 
epoxy had cured. 

Removed from the funnels, each sample was mounted onto one-half 
of its holder at each end (Figure 3.3) with epoxy. Before the epoxy cured, 
the mating half of the holder was bolted into place. This ensured that no 
slippage would occur within the holder. The gauge length used was 10 
inches. | 

To determine if interfiber friction, after a fiber failure occurs is 


significant some of the samples were lubricated with silicone oil. 


B. EQUIPMENT AND INSTRUMENT ATION 

The attributes of bundle testing are expedience in implementation 
and a substantial accumulatioin of information equivalent to a large 
number of single filament tests. 

To achieve the objectives of the experiment it was necessary to be 


able to apply a constant stress to the test samples. Figure 3.4, depicts a 
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Figure 3.1. Alcohol rinse setup. 
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Figure 3.2. Copper tube designed to produce laminar flow. 
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Figure 3.3. Holders shown in jig used to prepare samples. Gage length is 10 inches. 
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Figure 3.4. Block diagram of experimental setup. 
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block diagram of the instrumentation used. Its main feature is the IMC 
‘digital process controller. This processor (Figure 3.5) through automatic 
control of the movement of the Instron crosshead, using applied load as the 
control variable, creates a constant stress for a bundle with n variable 
number of fibers in time. 


For a viscoelastic material the governing intergral equation is 


O= f T E(t - tg) (de/de) dt (14) 


If E(t) were known, one could find de/dt and solve this equation for the 
desired stress o. But the problem was, that this equation could not be 
solved in real time. Additionally it may not even be linear. The solution, 
was to let the material solve this equation, whether it was linear or not 
by using a specially prepared control sample. 

The control sample was made from the same spool of graphite that 
the test samples were taken, and completely embedded in an epoxy resin. 
The reasoning behind this is the following: If a bare bundle of the same 
material were used as the controller, the compliance J(t,n) which is a 
function of time and the number of fibers would change. However, with the 
matrix a broken fiber is not entirely lost because the load is transferred 
Via the matrix binder and only a small ineffective length is lost. As a 
result J(t,n) will not be affected because the number of fibers remains the 


same. 
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Figure 3.5. IMC digital processor /controller. 
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Figure 3.6. Moving crosshead with load transducer. LYDT coil assembly provides 
displacement information to the {MC digital processor. 
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There are provisions in the Instron crosshead for nine samples. The 
center position is for the control sample. At this position (Figure 3.6) the 
control sample is attached to an Interface SSM-100 load cell which is 
attached to a 100 pound load transducer. The output of the load transducer 
is fed into the IMC digital processor. By setting the desired maximum and 
minimum load on the face of the digital processor and activating the 
automatic mode, the crosshead will move under the control of the 
processor to the desired maximum setting and continually adjust to correct 
for the effects of creep, within the-range set. 

Because the crosshead is rigid, all the samples of the same length to 
the left and right of the control sample experiences the same displacement 
and the same strain. Thus each fiber of the same length in the sample are 


subjected to the same strain and therefore, the same stress. 


o DATA ACQUSITION 

Data was recorded during the experiment using the Hewlett-Packard 
Interface Loop communications circuit (Figure 3.4). The HP-44468A Data 
Acqusition Control Pac provided the basic data logger software program. 
This program (Appendix A) was modified to permit the recording of time in 
hundredths of seconds, and edited to reduce the number of data registers 
used. It would have been desirable to have the HP-41CxX caculator 
manipulate the data while sampling, by having a comparison made between 
consecutive readings of the same channel, to account for changes in voltage 
as aresult of noise in the circuits. This was not feasible due to the design 
of the basic program and availability of data registers. 
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D. TESTING PROCEDURE 

Prior to the beginning of the test the 100 pound load transducer was 
calibrated in accordance with the manufacturers instructions. The load 
cells were checked for linearity and repeatability in voltage output, and a 
reading was recorded, for a weight of 2 kilograms for each load cell. 

For the first two experiments, three samples, one of which was the 
control sample and one of which was lubricated with silicone were place 
into the Instron, in positions 6, 7, and 8 (Figure 3.7). The crosshead was 
raised just enough to remove the slack in the three samples. By the use of 
a specially designed differential screw mechanism (Figure 3.8) each sample 
was individually loaded to 2 kilograms. With each sample at the same 
load, the equal length requirement of the experiment was satisfied. The 
crosshead was then lowered to zero load and the samples rested for a 
period of one hour to erase any history of stress before beginning the test. 

To start the experiment the IMC digital processors upper and lower 
limits were set to the prescribed load for that particular run. The 
crosshead control setting was adjusted to load the samples at the 
minimum rate possible. At the same time the loading was initiated, the 


data logger was activated. 
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Figure 3.7. Test samples loaded in the Instron. 
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Figure 3.8. Differential screw mechanism is designed to give a displacement 
of 0.0045 inch per turn. The upper portion is threaded 10-32, the lower 
portion is threaded 1/4-28 and both are right handed. 
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IV. RESULTS 


The experimental setup and procedure developed in this thesis 
produced acceptable initial results. However, the trial tests conducted 
have suggested that refinements in certain areas are necessary before 
actual testing begins. 

The placing of the samples into the Instron machine, and loading each 
individual sample to the same load was straightforward. The differential 
screw mechanism functioned as designed, providing adequate resolution and 
adjustment range for the graphite samples tested. However, a differential 
screw of different thread ratios may be required for samples with a 
drastically different compliance, as a consequence of different i iber type 
or different gauge length. 

The data acquisition rate was inadequate during the loading of the 
test samples at the start of the test. Installed in the HP-3421A data 
acquisition control unit were two 10 channel multiplexer boards, each with 
mechanical switching. The rate at which the unit switched from one 
channel to the next was fixed and relatively slow. During the initial 
loading sequence, it was important to obtain as many data points as 
possible in order to detect fiber breakage. In an attempt to compensate for 
the slow data acqusition rate the slowest possible setting on the Instron 
machine was used. 

During the first test, the load was applied at 18 grams per second, 


which appeared to be too fast for the rate at which the data was being 


ZS. 


recorded, but provided a reasonable contiuous load increase on the 
specimens (Figure B.1). The effect of inter-fiber friction is noted, by 
comparing the loading curves for the lubricated and the nonlubricated 
samples between the time interval of 4 to 7 minutes. In the dry bundle, as 
the fibers failed and became entangled they produced local load 
concentrations which caused premature fiber breakage. 

For the second test, the load was applied at a rate of 13 grams per 
second in an attempt to aquire more data points during the loading phase of 
the experiment. However, this was too slow, causing the drive motor of 
the Instron to stall as the load on the samples increased. The discontinuity 
In the load curves shown in Figure B.2a, for the first ten minute period 
resulted from the necessary readjustment of the loading rate on the Instron 
control panel in order to get the motor turning again. 

The IMC digital processor performed as expected. It maintained the 
constant load required on the test samples as evidenced by the bundle 
control curves shown in Figure B.1, for time greater than 9 minutes, and in 
Figures B.2a and B.2b, for time greater than 12 minutes. 

The effects of inter-fiber friction during creep rupture can be 
observed in Figures B.3 and B.4. The load supported by each individual 
bundle was normalized by their respective maximum loads attained at the 
peak of the loading ramp. As fibers fail as aresult of stress rupture, which 


is the stress per fiber times the number of fibers, the total load decreases. 


The normalized curve of F/R as represents the percentage of failed fibers 


In time. In fact, it Is the reliability function of the fiber. The silicone 


a0 


lubricated bundle shows a continuous decrease in supported load suggesting 
that the fibers failed sequentially, as expected. On the other hand, the dry 
bundle shows large discontinuities in it’s curve suggesting that inter-fiber 
friction within the bundle, caused by entanglement among the failed fibers 
caused high local stress concentrations leading to an accelerated increase 
In the number of fiber failures. 

During the constant load creep-rupture phase of the experiment, the 
time dependence of creep was physically observed by the operation of the 
IMC controller, causing the crosshead to move up in direction increasing 
displacement. The data recorded to support this observation was 
considered unreliable, in that it showed that the crosshead moved up and 
down. A polarity check of the LVDT and power supply used was conducted 
with no faults found. Had the excitation voltage been recorded, an 
explanation may have been possible. 

From the results of the trial experiments, some insight into the time 
dependent failure of graphite is possible. There are three mechanistic 
views with regards to the time dependent life of graphite-epoxy. One 
hypothesized mechanism is that graphite fiber is not time dependent, it has 
basically infinite life. What limits the life of the composite, is that weak 
fibers are broken during loading. The composite is sustained from failure 
by matrix load sharing which transfers the load to neighboring fibers. For 
a viscoelastic matrix such as epoxy, the spatial dimension required to 
transfer a given load, increases with time, thereby increasing the 


ineffective length, exposing additional weaker fiber sites to stress 
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concetrations which in turn cause additional fiber failures. Therefore, the 
composite failure mechanism is by Overload as a consequence of matrix 
creep. 

Another hypothesis is that the fiber filament itself has time 
dependent strength. This time dependent strength may be caused by flaw 
growth within the fiber as a result of macro-viscoelastic creep, cumulated 
from micro-slip among graphite slip planes, or by random bond breakage 
activated by stress or temperature. In either case, the failure mechanism 
of the composite is by time-dependent flaw growth of the fiber. 

Finally, the last hypothesis is that composite life is caused by both 
the viscoelastic deformation of the matrix and the time dependent flaw 
growth of the fiber. 

The preliminary results of the trial experiments may provide 
definitive evidence towards the conformation of the correct failure 
mechanisms. The life test results of the graphite bundles are presented in 
Weibull Coordinates in Figures B.S and B.6. Under such representation, a 
single mechanism of failure for a weakest-link configuration would 
manifest itself in a straight line. The test sample subjected to the lower 
stress level is represented in Figure B.S. This figure shows an essentially 
linear increase in failure probability as a function of time. This suggests 
that at lower loads the function is unimodal, and that fiber life is time 
dependent. Whereas, for the test sample subjected to the higher stress 
level (Figure B.6), there are distinct discontinuities in the curve’s slope 
occurring at 3 and 4.5 on the time scale, indicating that at higher !oads the 


failure density function is multimodal. This would tend to imply not 
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only that graphite life is not infinite, but that flaw growth mechanisms by 
stress activation and flaw growth mechanisms by thermal (or time) 
activation are different. Further investigation confirming these 
observations would be important for the characterization of graphite 
fibers for use in high performance applications and to provide input toward 


the improvement of fiber manufacturing technologies. 
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V. CONCLUSIONS AND RECOMMENDA TIONS 


A. CONCLUSIONS 

The experimental procedure presented proved to be a suitable method 
for lifetime testing of graphite bundles under constant load. The 
preliminary results of the trial experiments suggest that graphite life is 
time dependent and that the failure mechanism of flaw growth produced by 
temperature or time is distinct from that which is produced by stress. 
Through further refinement, the experimental technique developed offers 


promise of increased understanding of the failure processes. 


B. RECOMMENDATIONS 
The following are given as recommendations to improve the accuracy 
of the data. 
l. Data Acquisition: \nsert the HP-85 data logging system into the 
HP-IL loop for the initial loading of the bundles. This unit can record 
data at a much faster rate. 
2. Sample Freparation: Lubricate the test samples with silicone oll 
to eliminate the catastrophic effects of inter-fiber friction. 
3. Frocedure: Load the samples at a minimum of 18 grams per 


second. This rate will prevent the drive motor from stalling. 
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FigureA.!. Modified HP44468A Data Logger Routine (cont'd) 
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Figure A.2. Modified HP44468A Output Format Routine 
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Figure B.1. Test 1: Load-Time curves. Maximum constant 
load on the control sample was 12.575 kilograms. 
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Figure B.2a. Test 2: Load-Time curves. Maximum constant load on the control 
sample was 16.003 kilograms. 
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Figure B.2b. Continuation of Test 2: Load-Time curves. 
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Figure B.3. Test !: Normalized load curves. Maximum load 
achieved by the bundle with silicone was 8.430 kilograms. 
Maximum load achieved by the dry bundle was 8.900 kilograms. 
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Figure B.4. Test 2: Normallized load curves. Maximum load 
achieved by the bundle with silicone was 11.572 kilograms. 
Maximum load achieved by the dry bundle was 11.329 kilograms. 
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Figure B.S. Test 1: Lifetime distribution curve. 
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Figure B.6. Test 2: Lifetime distribution curve. 
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